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ABSTRACT 

The naval importance of knowledge of short-period variations in 
ocean thermal structure is exemplified by the effect of internal waves 
upon underwater sound propagation patterns. This thesis describes an 
investigation undertaken to examine some of the characteristics of 
short-period (.071-.01 ae eet variations of thermal structure in 
the vicinity of Monterey Bay. Time series bathythermograph obser- 
vations at 7.0 minutes'interval and to ten hours duration were col- 
lected at a position in 65 fathoms near the shoulder of the Monterey 
submarine canyon. The depths of one-half degree Fahrenheit isotherms 
were analyzed for power spectral distribution, and comparison of the 
individual time-series'spectra was made. Variance of depth of selected 
isotherms over 140-minute intervals was compared with tidal phase, 
showing maximums of variation occurring immediately after tidal extremes. 
A comparison was drawn between spectral analysis of data obtained from a 
thermistor temperature sensor suspended in the water beneath a pier in 
Monterey Harbor. Internal wave-like characteristics are recognizable, 


though conclusive evidence of the existence of internal waves cannot be 


shown from the data collected. 


alga 








TABLE OF CONTENTS 


Section Title 
Ll. Introduction 
2. Program Objectives and Description 


a. 


le 


Objectives 

Observation Site Selection 

(1) Offshore Location 

(2) Nearshore Location 
Equipment and Procedures 

(1) Stationing 

(2) Bathythermograph Employment 
(3) Miscellaneous 

(4) Nearshore Observations 


(5) Data Processing 


Description of Data 


ad. 


|e} 


Cc. 


d. 


Data Collected 16 March 1965 
Data Collected 26 March 1965 
Data Collected 2 April 1965 


Data Collected 6 April 1965 


Discussion of Analysis and Results 


a. 


[Be 


Cc. 


d. 


Power Spectrum Analysis 
Modes of Oscillation 
Tidal Phase Relationship with Isotherm Variance 


Comparison of Offshore and Nearshore Spectra 


Summary of Conclusions 


iii 


Page 


10 
2 
iS 
19 
21 
Pag 
AJ 
ou 
32 
34 


38 








Figure 


10 
Wa 
Wa 
13 


14 


LIST OF ILLUSTRATIONS 


Chart of Vicinity of Observations 

Plot of Raw Data; 16 March 1965 

Typical BT Traces; 16 March; 2 April; 6 April 
Plot of Raw Data 26 March 1965 

Typical BT Traces; 26 March 

Plot of Raw Data; 2 April 1965 

Plot of Raw Data; 6 April 1965 

Power Spectra; 16 March 1965 

Power Spectra; 26 March 1965 

Power Spectra; 2 April 1965 

Power Spectra; 6 April 1965 

Power Spectra by Summation of Bands 
Isotherm Variance Related to Tidal Phase 


Power Spectra (Nearshore and Offshore) 


Lv 


Page 


13 
14 
16 
ly 
20 
22 
25 
24 
jis, 
26 
30 
33 


35 








1. Introduction 

Much effort has been directed toward the description and prediction 
of the long-period variations in ocean thermal structure for the purpose 
of estimating variations in underwater sound propagation velocity. 
Methods of predicting characteristics of the thermocline, to which many 
sonars are sensitive, have been developed and tested (Tabata [13] : 
Laevastu [12] ). The usefulness of such techniques is limited, however, 
in that they provide, at best, an estimate of mean characteristics about 
which some variation can be described. These authors suggest that a 
portion of the variation can be explained as the effect of internal waves. 

Lee [8] has shown a dramatic example of the calculated sound re- 
fraction, and the resultant varying-intensity sound energy field, due to 
the presence of idealized high-frequency internal waves representative 
of those observed in nearshore areas off San Diego. It is clear, from 
his example, with its areas of marked intensification and attenuation, 
that the effect of internal waves is more complex than a uniform change 
or variation in thermocline depth. Additionally, it is apparent that 
the single bathythermogram cannot truly describe the thermal structure 
of a vicinity, in that the nature of variability is not indicated. Thus, 
there seems ample justification for the pursuit of more knowledge of the 
short-period variations in thermal structure, ultimately, of course, for 
the evaluation of their effect upon sonar detection range. This project 
was undertaken to provide a description of some such thermal structure 


variations. 








2. Program Objectives and Description 


a. Objectives 
The project was planned as an investigation of short-period 
oscillations of internal temperature structure in the vicinity of 
Monterey Bay, hopefully to show the existence of internal waves and to 
describe their characteristics. To this end, simultaneous observations 
of subsurface temperature variation were planned at two sites, one off- 
shore using a bathythermograph, the other nearshore, employing two ex- 
isting thermistor temperature sensors with recorders installed on Wharf II 
in Monterey Harbor. The nearshore site was not ideally located for 
observation of internal waves progressing from offshore; however, the 
spectral distribution of the variations observable from the pier would 
be of interest. It was intended that several time-series would be col- 
lected, each of duration (approximately) of the tidal period. The 
several series were to be collected during tides of various magnitudes 
so that some evidence of the effect of tidal current variations upon 
internal wave generation might be noticed. 
Observation of a deep scattering layer, if any appeared, was planned 
to determine the effect, if any, of short period thermal variations 
upon its depth. 
b. Observation site selection 
(1) Offshore location 
The offshore location (figure 1) was selected, at a position 
approximately 36°40'N, 122°00'W (292T/3% miles from the Pt. Pinos buoy), 
after consideration of several factors, including evidence of both tem- 


poral and horizontal variation of vertical thermal structure. 
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The considerations follow: 

(a) It was near the position selected that Skogsberg (1] 
was unable to duplicate temperature measurements at depth within periods 
as short as one hour. 

(b) In July 1962 the USS Marysville (EPCER-857) towed a 
thermistor chain in the vicinity of the location selected, and recorded 
(record courtesy of Dr. G. H. Jung) variations in depth of isotherms 
ranging from 30 ft. in relatively steep gradients, to over 100 ft. where 
gradients were not severe. "Apparent wavelengths' ranged from 1/8 mile 
to three or more miles respectively. 

(c) In that several authors (Davis and Patterson [10] , 

Lee [9] ) have suggested in theory or shown through observations and 
model studies that internal waves may be formed as tidal currents flow 
over bottom discontinuities, the site was positioned on the shoulder of 
the Monterey submarine canyon. The water depth is approximately 65 
fathoms, and the bottom slopes steeply seaward, but gradually shoreward. 
The bottom contours are somewhat convex seaward, so may provide con- 
vergent refraction of internal waves moving shoreward to the site. The 
bottom is sufficiently level, at the site, so that likelihood of anchor 
movement on slopes is reduced. The water is sufficiently shallow so that 
the entire water column was subject to observation with reasonably high 
frequency of BT drops. 

(2) Nearshore Location 

The nedpanene observation site was determined by the location 
of existing thermistor temperature sensors located beneath the Wharfinger's 
Office on Wharf II in Monterey Bay. The progress of internal waves from 


the seaward observation position involves refraction through a minimum 
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of 50°-60°, (see figure 1) and travel-distance a minimum of seven miles. 
c. Equipment and Procedures 

The offshore observations were conducted from the U. S. Naval 
Postgraduate School Oceanographic Research Vessel (63 ft.), using a 
450 ft. BI with gold-wash BT slides. The boat navigated to within 4% mile 
of the selected site at 65 fathoms where a taut-line moored marker buoy 
was then fixed. Position of the boat was maintained, using engines, by 
visual reference to the marker buoy. Several visual and radar fixes 
during each observation period assured no drift of the marker. 

(1) Stationing 

The taut-line moored marker was employed since no anchoring 

gear available was satisfactory for anchoring the boat. The moored 
float rig was designed, with reference to Isaacs, et al [4] , to mark 
the anchor position to within 100-ft. radius. It consisted of a scrap 
iron anchor in an "H" form weighing approximately 120 pounds in air, 
1/8 inch diameter (No. 4) braided nylon mooring line, and a float con- 
structed of a commercially available 13-gallon capacity plastic carboy. 
The marker buoy was placed on station immediately prior to each obser- 
vation period. When fathometer-indicated depth reached 65 fathoms 
during the approach to station, the boat was stopped and the marker- 
buoy anchor lowered to the bottom using the nylon line over a gypsy-head 
on the winch. The boat was maneuvered to the anchor position using 
"wire-angle'" as a guide, and the nylon line hauled in by hand until more 
than 30 pounds (estimated) of tension had been established. The float 
was then tied on and allowed to drift free of the boat, marking the 


anchor position. 








An alternate method of lowering the anchor was used for 
safety, particularly, in adverse seas. A one-inch manila line was 
used to lower the anchor to the bottom, while the nylon line attached 
to the anchor was paid out slack to avoid damage and strain. Once 
bottomed, the anchor was raised approximately ten feet from the bottom, 
the nylon line was tensioned as before, by hand, and the float attached. 
The anchor was then lowered once again with the marker float in the 
water. When the marker float was observed to be operating satisfactorily, 
its draft a rough indication of sufficient tension in the mooring line, 
the manila lifting line, with ample slack, was attached to a second 
supporting float consisting of three plastic one-gallon Clorox bottles, 
and allowed to drift free for use later during recovery. The manila 
lifting line, thus supported, provided little if any lifting force on 
the anchor and insufficient dragging force to move the anchor. 

The No. 4 nylon provided sufficient margin of strength for 
the operation, yet was sufficiently elastic through the range of loads 
imposed by surface waves, wind and current, to maintain nearly constant 
tension. It cannot endure much abrasion damage, however, and care was 
taken in handling. 

Station-keeping required that the helm be manned constantly, 
and difficulty varied considerably depending upon wind and swell relative 
direction and intensity. Long period waves (swell) to eight feet, and 
winds to approximately 18 knots were limiting maximums. [In general, 
station was maintained to within 100 yards, usually to within 100 ft. 
of the marker. When swell and wind direction were nearly constant during 


an observation period (the usual situation) the direction of error from 








on-station varied little. The buoy was kept "close on the bow", and the 
bow headed constantly into the wind, which supplied the major drifting 
force. 
(2) Bathythermograph Employment 

Once position had been established, BT drops were commenced 
at seven-minute intervals, and continued until the practical limit of 
endurance, equipment breakdown, or adverse weather conditions. The 
interval was established after determining through brief experimentation 
that good BT traces were obtained when a raising/lowering rate of 50 
meters/minute was used. At this rate, neither "double-traces" due to 
lag of the BT mechanism during stationary dropping nor "stepping" due 
to static friction effects were noted. The gold-wash BT slides exhibited 
a tendency to flake, occasionally, necessitating cleaning of the stylus 
after each drop to remove residue which might produce a wide trace on 
subsequent slides. A controlled hydrographic winch with '"meter-wheel" 
was used to provide the selected lowering/raising rate. A 150-pound 
weight was suspended ten feet beneath the BT to provide stability and 
to aid in making drops vertical. The seven-minute data interval allowed 
sufficient time between drops for changing slides, cleaning stylus, and 
allowing at least 30 seconds of soaking to stabilize the BT at the sur- 
face prior to each drop. 

(3) Miscellaneous 

Occasional observations of "bucket temperature", wind, wave 
conditions, and cloud cover were recorded during each period. Several 
measurements of salinity, conductivity, and temperature were attempted 


at 50 ft. depth intervals to 300 ft. using an Industrial Instruments, Inc. 





RS5-3, Portable (induction) Salinometer. However, operation was not 
entirely satisfactory, and depth of sensor impossible of accurate 
determination due to streaming in relative currents; hence, data ob- 
tained were not reliable. 
(4) Nearshore Observations 

Temperature variation with time was measured by two therm- 
istor temperature sensors, with appropriate amplifiers, and two Varian 
pen recorders operating with paper speed of one inch/hour. Full-scale 
deflection of the record was 5°C on five-inch recording paper. The 
sensors were suspended from Wharf II in Monterey Harbor in 28 feet of 
water at nominal depths of eight feet and 15 feet, which varied as water 
level did. Maintenance difficulties prevailed in all but one offshore 
observation period when a satisfactory record was obtained during a 
period of 18 hours and 40 minutes, but lagging the offshore observations 
by about nine hours. 

(5) Data processing 

The BT slides were viewed through a standard handheld 
magnifying viewer and depths of occurrence of each 4%F noted. When 
double-tracing occurred, the depth selected was the arithmetic mean of 
the two depths of crossing. A wide section in a portion of the trace 
was treated as incipient double trace. When a wide trace appeared 
throughout the BI slide record, it was generally interpreted as due to 
debris on the BT stylus, and the center of the trace was estimated as 
the correct indication. These procedures provided data generally of 
very good continuity when integrated with the rest of the series. (The 


author personally handled all BT drops and extracted all data, so 
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consistency of procedure and accuracy of data transcription is believed 
good). It should be noted that accuracy of interpretation varied as the 
temperature gradient, less depth accuracy possible when gradient was 
"weak'"', In general, independent reinterpretation of selected BT traces 
produced data within three feet of the original estimate. 

All data were then plotted as isotherm depth versus 
time to check apparent continuity and reasonability of the time series 
data accumulated. Thus checked, the data of each isotherm were pro- 
cessed for power spectral estimates on a Control Data 1604 computer 
using the BIMD-35 (Power Spectral Analysis) Program, one of a series of 
programs compiled by the University of California, Los Angeles, Division 
of Biostatistics, Department of Public Medicine and Health, School of 
Medicine. The program was employed in several ways, i.e., using 
various prewhitening constants, removal of linear trend and using 
various numbers of lags in order to gain experience. The analysis used 
to gain the results presented, however, is without prewhitening or 
detrending, and employs arbitrarily selected numbers of lags. 

Data? taken at Wharf II, the nearshore site, were 
not plotted in raw form, since the recorder trace offered the visual 
information necessary for estimating the validity of the data. Tempera- 
ture was read from the trace at seven-minute intervals and placed on cards 
for computer processing in the same manner as the offshore data for com- 
parison of spectral distribution with that estimated from the appropriate 


data collected offshore. 


1tnstrumentation on Wharf II supplied under research grant to U. S. 
Naval Postgraduate School by the Office of Naval Research. 





3. Description of data 

BT data were collected on seven occasions at the designated site. 
Four of the sequences obtained were long enough to warrant spectral 
analysis. The longest of the three rejected extended only three hours. 
BT slides and numerical data extracted from all of the sequences are in 
the files of the Department of Meteorology and Oceanography, U. S. Naval 
Postgraduate School. This section describes the four significant 
sequences and the conditions surrounding their collection. The data 
are plotted in raw form as 4F isotherm depths versus time of observation. 

Some of the series indicate isotherm oscillations of relatively 
large magnitude where temperature gradients are weak. It should be 
noted that error of interpretation and of the instrument is magnified 
in this sort of presentation when the gradient is weak. The BT trace 
may contact the temperature index line for 50 feet or more of depth and 
estimates of the exact depth of crossing become poor. 

There are some series of oscillations which are described by only 
one data point per extreme. One might suspect instrumentation noise 
since the wave period is twice the observation interval. It can be noted 
in some of these incidents, however, that the isotherms at other depths 
are progressing smoothly in longer period oscillations, so support is 
lent to the credence of data at all levels. 

To evaluate the effect of positioning error upon the significance 
of data, assuming that it represents internal waves, some wave lengths 
were calculated. A constant density gradient was assumed since a marked 


thermocline does not exist in most instances observed. Salinity and 


10 





temperature data collected allowed calculation of an approximate density 


gradient of -3.6 x 107° g/cm>/m. The wavelength for this condition 


is given by the formulas: 








Ch = C™ = h n2- o 2 
Kn AAW 
= an 
L = 2 ‘VW Ge 
on 
where: 


n = mode, an integer 


Cc, = phase speed 
GO =2% = angular frequency 


iD 
h = water depth 


L = wave length 


N =/- g @p = Vaisalla's frequency 


if Oz 


(3208) 


(3507) 


The wave length of the shortest waves observable (.071 cpm) from the 


data is greater than 540 meters, assuming a mode-1 oscillation (a maximum 


of vertical velocity occurring at only one level in the water colum.) 


The variance of position was generally less than one quarter the wave- 


length of the shortest observable feature. 
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a. Data collected 16 March 1965 


The following describes conditions surrounding the collection 


of the raw data plotted on figure 2: 


Position 


Depth of bottom 

Time between BT drops 

Times of observation 

Length of observation period 
Number of data 

Wind condition 

Estimated swell 


Sky conditions 


Tides 


Comments: 


26040.15'N 
121°59.80'W 

59 fathoms 

Seven minutes 
0845-1606 PST 

7 hours, 21 minutes 
64 points 

2102/6" kt. 

270°/3-4 ft. 


Overcast with occasional brief 
breaks (ten minutes maximum). 


Low 0346 PST/1.5 ft. 

(Range +3.9 ft.) 
High 0915 PST/5.4 ft. 

(Range -6.0 ft.) 
Low 1603 PST/-0.6 ft. 

(Range +5.5 ft.) 
High 2221 PST/4.9 ft. 


A typical BT trace, figure 3, shows no marked thermocline. Some double 


tracing of the BT occurred, especially near mid-depth and during the 


latter 2% hours of the sequence. 


The lack of sunshine during the day 


is reflected in little warming indicated, the near surface isotherm 


trending downward only slightly during the period. 


Density calculations from salinity data collected at several depths 


between the surface and 250 feet, indicates density gradient of approxi- 


mately -3.6 x lore g/cm3/meter. This value is an estimate only, since 


the data are not accurate. 
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b. Data collected 26 March 1965 


The following describes conditions surrounding collection of 


the data plotted on figure 4: 


Position 


Depth of bottom 

Time interval of BT 

Times of observation period 
Length of observation period 
Number of data 

Wind condition 


Estimated swell 


Sky conditions 


Tides 


Comments: 


36940 .20'N 

122900 .05'W 

65 fathoms 

seven minutes 
0744-1753 PST 

10 hours, 9 minutes 
88 points 


Calm 


Initially 280°/3 ft., 6 second period 
1600 PST 280°/2-3 f£t., 10-13 second 


period 


Overcast initially, with occasional 
light showers. Clouds thin and broken 
during early afternoon, then thin 


overcast after 1530 PST 
High 0507 PST/4.0 ft. 
Low 1231 PST/0O.2 ft. 


High 1927 PST/3.8 ft. 


(Range -3.8 ft.) 


(Range +3.6 ft.) 


The thermocline which was observed early in the period at about 75 ft. 


persisted through the morning, appeared to separate into two "steps" by 


early afternoon, and dissipated entirely by the end of the series. Many 


double traces of the BI were noted, some in the region of the thermocline, 


suggesting oscillations with periods less than sampling interval. A 


typical BT trace is shown on figure 5. 
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Data were collected at Wharf II during this period, but the useful 
record obtained extended only from 1600 PST, 26 March 1965, to 1040 PST, 
27 March 1965. The temperature record obtained from the 8-feet deep 
sensor appeared reasonable, with a temperature range of 2°C. The 15-feet 
deep sensor is believed to have been operating, but the record indicates 
only slow temperature drift through a range less than 0.25C. The pier 
record is not shown, but its power spectrum has been calculated (8 ft. 


sensor only) and is discussed later. 
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c. Data collected 2 April 1965 


The following describes conditions surrounding collection of 


the data plotted on figure 6: 


Position 


Depth of bottom 

Time interval of BT 

Time of observation period 
Length of observation period 
Number of data 

Wind condition 


Estimated swell 


Sky conditions 


Tides 


Comments: 


36°40.20'N 
122°00.05'W 

63 fathoms 

Seven minutes 
0744-1828 PST 

10 hours, 44 minutes 
93 points 

Calm (NW) 


Initially 3 ft/10 sec/280°T 
1030 PST 6 ft/10 sec/280°T 
Decreased slowly thereafter 


Occasional brief sun through breaks in 
thin overcast. Several short periods 

of drizzle until early afternoon when 

clouds became scattered (.5 cover) 


Low 0438 PST/0.5 ft. 

(Range 3.8) 
High 1042 PST/4.3 ft. 

(Range -3.3) 
Low 1644 PST/1.0 ft. 

(Range 3.8) 
High 2243 PST/4.8 ft. 


Transient thermoclines appeared and disappeared during period. Double 


tracing was common, especially during the last several hours of the 


sequence. Minor inversions and inversion tendencies appeared, often 


accompanied by double traces. 


figure 3. 


A typical BT trace is plotted on 
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d. Data collected 6 April 1965 


The following describes conditions surrounding collection of 


the data plotted on figure 7: 


Position 


Depth of bottom 

Time interval of BT 

Time of observation period 
Length of observation period 
Number of data 


Wind condition 


Estimated swell 


Sky conditions 


Tides 


Comments: 


36°40 .23'N 
121°59.50'W 


66 fathoms 

7 minutes 

0713-1715 PST 

10 hours, 2 minutes 

87 points 

Initially 280°T/3 kt. 
1100 PST - 285°T/8 kt. 
1300 PST - 290°T/15 kt. 


Initially 2 £t/280° 


0933 PST - occl. 5 ft. swell 
1200 PST ~ 1-2 ft. wind waves on swell 


1500 PST ~ sea 310T/7 ft. 


Overcast with occasional breaks during 
morning, thinning during afternoon 


High 0013 PST/1.9 ft. 
Low 0751 PST/-0.4 ft. 
High 1455 PST/3.4 ft. 


Low 1922 PST/2.8 ft. 


Many small transient features were observed on the traces. 


double traces occurred, they were not prevalent. 


(Range -5.3 ft.) 
(Range 3.8 ft.) 


(Range ~-.6 ft.) 


Though 


Surface warming during 


the afternoon is indicated by down trend of near surface isotherms. 


Typical BT trace is shown on figure 3. 
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4. Discussion of Analysis and Results 


The observed variation in depth of isotherms is treated in parts of 
this discussion as though due in significant part to the vertical move- 
ment of water in internal waves, during which water-particle temperature 
is conserved. It is recognized that other effects, not easily isolated, 
may be present also, as the advection of inhomogeneous water, convective 
and current-driven eddies, turbulence due to wind and wave action, 
radiative and molecular heat exchange, chemical reactions and sampling 
error. Suspicion that internal waves are responsible for a large 
portion of the temperature oscillation is enhanced by the apparent 
coherence in phase between isotherms through a major portion of the 
water column (see figures 2, 4, 6 & 7). It is more difficult to assign 
the cause of such activity, that is oscillations apparently random in 
period and amplitude yet coherent in large vertical extent, to other 
phenomena such as eddies. (If coherence were notable in all isotherms, 
then misplacement of the slide in either the bathythermograph or on 
the viewing grid could be responsible.) Regardless of which phenomenon 
it is that contributes most to the observed fluctuation, the spectral 
analysis can provide some insight regarding the variability of thermal 
structure as observed by the BT. 

a. Power Spectrum Analysis 

The results of the power spectrum analysis of isotherm oscil- 
lations are on figures 8 thru 1l. The spectra of continuous OF iso- 
therms from each observation period were plotted in an array ordered 
according to the median depths of isotherms, so a rough impression can 


be gained of the position in the water column which each spectrum 


se 





describes. The power spectrum estimates were computed employing 

various numbers of lags. (See Appendix I ). The number of lags 
employed for the calculation of the spectra presented is noted on 

each of the figures. The number of lags chosen exceeds the 54-10% of 
the record recommended by Blackman and Tukey (see Jenkins [7] ) as the 
appropriate truncation point. Lags used, rather, comprise 25% to in 

one case 38% of the length of the record. The resultant power spectrum, 
even with such a large number of lags, appears quite smooth and exhibits 
characteristics much like the spectra calculated with fewer lags. This 
is apparently a reasonable result since, as suggested by Jenkins, 
truncation points may be much larger than those suggested by Tukey if 
the autocovariance function converges to negligibly small values at 
greater lags. It does become small (less than 10% of zero lag auto- 
covariance) in several of the series analyzed with 35 lags. No attempt, 
however, has been made in this project to optimize the spectral analysis 
through selection of various lags as Jenkins proposes. More meaningful 
results could be described were such techniques employed. 

The presence of long-period oscillations presents a problem in the 
analysis. The data records are shorter than the period of low frequency 
harmonics which can be observed in the raw data plots. These long period 
contributors may be of tidal period or local inertial period (10.04 hrs.). 
Such an effect interferes with the estimation of the spectrum at other 
frequencies so that removal of long-period influence by “detrending" may 
be of some help. Detrending was employed experimentally; but, again, 
negative spectral estimates were sometimes obtained (see Appendix I ), so 


computations were made without detrending. A better solution would be 
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precise removal of the long period oscillation, using a best-fit 
sinusoid determined from a longer series of data with appropriate 
smoothing. 

A visual comparison of area under the power curves suggests that 
more energetic oscillations occurred during the period of largest tidal 
range. Maxima and minima of power apparent on the curves cannot, in 
most cases, be shown statistically significant with 90% confidence by 
the chi-square test (Panofsky [11] ). However, some of the power 
extremes are supported by the appearance of similar features in the 
power spectra at several adjacent levels in the array. By adding the 
area under all the isotherm spectrum curves of all four observation 
periods by 0.0l-cycle/minute bands, it is shown in figure 12 that the 
mean power spectrum is a smooth one with no frequencies apparently 
favored. The maxima indicated in the arrays of individual isotherm 
spectra, then, are either the result of sampling fluctuations or of 
internal wave spectral distributions peculiar to the conditions during 
the observation period, but none are permanent features of the mean 
smoothed power spectrum. 

Estimates of power spectra are, in Sell, associated with the 
assumption that the time series represents a stationary process. In 
this application, the power spectrum is used to describe short-period 
oscillations, which may be internal waves, suspected of formation in 
some way associated with tidal phenomena. The tides, however, do not 
represent a steady—state process; so the assumption cannot be satisfied, 
except approximately, for short segments of the tidal cycle; and, un- 


fortunately, the spectrum cannot be found accurately and with definition 
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from a segment of short duration. In addition to the effect of con- 
tinuous variation of the suspected exciting agency, the effect of in- 
ternal waves near the period of the tide tends to interfere with esti- 
mation of the power in the higher frequencies of interest. Long records 
permit better estimates of the mean power spectrum during the interval; 
but,if the record includes all stages of the tide, nothing can be learned 
about the relationship between high frequency energy and tidal stage. 
Perhaps with many shorter segments of data, carefully selected with 
regard to similar conditions of tidal phase, amplitude, and water 
physical structure, some more significant relationship could be de- 
termined by application of the ergodic theorem (see Davenport and Root [5] De. 
At any rate, optimizing the techniques of spectral analysis for ex- 
ploratory investigation of relationships between tidal activity and 
high frequency internal wave activity would require mathematical 
treatment beyond the scope of this project. 
b. Indications of Modes of Oscillation 

Examination of the isotherm oscillations (figures 2, 4, 6 & 7) 
reveals a suggestion of higher than first-order oscillations in some 
regions; but, since a spectrum of frequencies contributes to the oscil- 
lations observed in the raw data, clear evidence of higher modes is not 
present. A continuous distribution of density with depth is suggested 
by the data collected. Internal waves, therefore, are likely to exist 
in a variety of modes of oscillation (Defant [6] ). The first mode of 
oscillation is that in which vertical displacement,or vertical velocity, 
is in the same direction throughout the water column, with the maximum 


occurring at one depth. Second mode describes opposite directions of 
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motion in upper and lower portions with two maxima of amplitudes; third 
phase has three maxima, etc. 

The power spectra (figures 8 thru 11) may also be examined for 
suggestions of higher modes of oscillation. The spectral distribution 
represents an indication of the contribution to total variance by 
oscillations in particular frequency bands. A comparison of the mag- 
nitude of the power function in a particular frequency band at various 
levels in the water column, then, suggests levels at which the variance, 
or the amplitude of vertical isotherm movement, is greater. These may 
also be interpreted as levels of antinodes, indicating the mode of 
oscillation present. There are some instances of maxima of the power 
function at more than one level in the arrays of power spectra, as 
near .018 cycles/minute on figure 11. This might imply energy in a 
higher mode; but again the evidence is not conclusive. Sampling error 
cannot confidently be discounted as the cause for apparent spectral 
maxima. Additionally, if the process is non-stationary, then oscil- 
lation in a particular mode of a particular frequency may be unlikely 
to persist throughout the observation period. 

c. Tidal Phase Relationship with Isotherm Variance 

A gross indication of internal wave activity relative to the 
phase of the surface tide wave is shown on figure 13. The variances of 
20-point segments of the data (about each 20-point mean) were calculated 
for two isotherms from each of the four observation periods. Isotherms 
chosen were those nearest depths of 100 feet and 200 feet. In most 
cases, the magnitude of variance appeared Peedlwet in the segments of 


data which represent the period immediately following the tidal extremes. 
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Some relationship between tide and high-frequency thermal fluctuations 
is suggested by this nearly consistent association of variance and tidal 
phase, though in some cases the maxima of variance are quite small. 
Knowledge of the relationship between tidal phase and tidal current in 
this region would allow more interesting speculation upon the causal 
effect of tidal currents. 
d. Comparison of Offshore and Nearshore Spectra 

The relating of offshore and nearshore data by coherence and 
cross-spectra was not accomplished due to computation difficulties. 
Processing of nearshore data consisted merely of the calculation of the 
power spectrum of temperature oscillations with time. The only period 
of nearshore data which corresponds with offshore observation lags the 
offshore period by 8 hours 16 minutes and extends for 18 hours 40 
minutes. 

The power spectrum of the nearshore data is shown on figure 14. 
The power spectrum of the offshore isotherm (51.0°F) representing nearly 
the same relative depth is also plotted for comparison. The smoothed 
shape of the spectra are similar, and certain similarities in detail 
may be noted too. Of all the isotherm spectra from offshore data the 
one shown is most like the nearshore spectrum in both smoothed shape and 
detail. The chi-square 90% confidence limits are also shown. (Other 
isotherms in the same series appear in part on figure 9). 

There is little possibility that advection phenomena could 
account for the appearance of similar characteristics offshore and 
nearshore. An investigation of currents in the vicinity of the near- 


shore sensor (Breidenstein and Thomas [3] ) suggests that they are 
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primarily due to tides with secondary surface wind effects. It was 
mentioned earlier that the temperature at 15 ft. depth nearshore varied 
little, while the temperature at 8 ft. depth varied considerably; this 
suggests internal waves rather than advection phenomena. 

Calculation of the phase speed of a high-velocity internal 
wave, using the data from the offshore site and equation (3.01), yields 
0.702 meters/second (1.36 kt.). If progressive internal waves originate 
offshore, then it is possible, accounting for the reduction of velocity 
assuming linear shoaling, that some of the energy observed offshore is 
observable in the nearshore data. Assuming that group velocity equals 
phase velocity, the first of the highest velocity waves discernable in 
the data (those of Nyquist frequency) would arrive, at the earliest, 
5.2 hours after the beginning of the nearshore data series. Energy 
observed offshore may be observable in no more than 13.5 hours of the 
nearshore data. 

The apparent similarity of the spectra suggests that they 
represent the same type of phenomena. The same source may produce all 
the perturbations observed, or a similar mechanism may be active but 
in different source areas. Calculation of coherence between data from 
both locations would aid in determining whether the same source produced 
them. Comparison of the spectra does not indicate dispersion of energy 
due to the frequency dependence of internal wave velocity. One might 
anticipate some filtering of wave frequencies due to dispersion in a 
single train of waves progressing shoreward from the offshore site, the 


predominant frequencies depending upon the lag of observation. The 
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absence of filtering may suggest that internal waves are generated con- 
tinuously, or that a similar generating mechanism produces the phenomenon 
near both positions. 

The justification for accepting or ignoring a suggestion that 
offshore phenomena may be propagated or advected to a nearshore position 
could be established. Determination of the nature of the phenomena 
would be necessary. If offshore temperature data were collected simul- 
taneously at three locations (triangular), then direction and speed of 
propagation of the fluctuations could be calculated. Mean current 
measurements would allow discrimination between advective or internal 
wave phenomena. Such discrimination might be made by noting whether 
the disturbance moved with the velocity of the current, or possibly 
relative to the current with velocity consistent with that of internal 
waves. Even more detailed current observations might indicate hori- 
zontal wave particle velocities (periodic currents) coherent with the 
vertical motions; this would be most significant evidence of the 
existence of internal waves. Such detailed current measurements would 


be difficult to make and analyse in such a spectrum of perturbations. 
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5. Conclusions 

a. The spectral distribution offshore and nearshore appear similar 
in smoothed form in the one available instance of comparison. 

b. A relationship between tide phase and high frequency oscillations 
is suggested by the increase of variance of isotherm depth immediately 
following the tidal extremes. 

c. Spectral peaks, if real, appear to be transitory in nature both 
offshore and nearshore. No characteristic detail persists through all 
ovservations. 

d. Conclusive evidence of the existence of internal waves was not 
shown and the progression of isotherm perturbations with the character- 
istics of internal waves could not be determined from the data of the 
two isolated observation positions. 

e. Integrated spectral density appears greater in the presence of 
larger amplitude tidal fluctuations. No quantitative estimate is offered. 
f. The application of power spectral analysis appears a promising 

tool in the investigation of internal wave generation relationships, 
but is not optimized in this project. 

g. Accurate-data collection remains an unpleasant obstacle in the 


pursuit of oceanographic knowledge. 
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APPENDIX I 

Initial test computations with large numbers of lags (as many as 
50% of the number of data) were performed, resulting in some unexplained 
negative power estimations apparently due to instability in the numerical 
procedures. Similar difficulty was encountered when attempting to cal- 
culate the cross-spectra and coherence of offshore and nearshore data, 
so results of this sort are not presented. Professor R. Read, U. S. 
Naval Postgraduate School Mathematics Department, believes that this 
unexpected result may be due to the fact that the method used in the 
BIMD-35 program for estimating the autocovariance function does not 


force the estimate to be a positive definite function. f(p) is a 


pesitive definite function if for all choices of x], x9,..... 3. 2+Xp> 
nw” 
) Lx £G-mm, 2 0° 
jel het 


The number of lags employed for calculation of results shown were lower 


than those yielding negative spectral estimates. 
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